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ABSTRACT 

Combat attrition can be modelled as either a deter- 
ministic or a stochastic process. In this thesis the fore- 
casts of combat outcomes generated by deterministic and 
stochastic Lanchaster-type models for combat between two 
homogeneous forces are compared. Using stochastic formula- 
tions, the probability of winning is studied and contrasted 
with a deterministic win criterion for some idealized combat 
situations. Force ratios required to assure victory in an 
attack are studied from this standpoint. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION 5 

A. REVIEW OF LANCHESTER THEORY OF COMBAT 5 

B. STUDY OBJECTIVES 6 

II. COMBAT BETWEEN TWO HOMOGENEOUS FORCES 8 

A. LINEAR LAW ATTRITION ON BOTH SIDES 8 

1. Assumptions 8 

2. Deterministic Model 9 

3. Stochastic Model 10 

B. SQUARE LAW ATTRITION ON BOTH SIDES H 

1. Assumptions H 

2. Deterministic Model H 

3. Stochastic Model 1^ 

C. THE MIXED CASE H 

1. Assumptions 1^ 

2. Deterministic Model 1^ 

3. Stochastic Model 

III. COMPUTATIONAL RESULTS 

A. LINEAR LAW 1^ 

B. SQUARE LAW 

C. MIXED LAW 

IV. SUGGESTED FUTURE WORK 

V. CONCLUSIONS 



3 



f 










COMPUTER PROGRAM 



38 



LIST OF REFERENCES 40 

INITIAL DISTRIBUTION LIST 41 

FORM DD 1473 42 



4 



4 

% , 




1 ‘ 







I. 



INTRODUCTION 



A. REVIEW OF LANCHESTER THEORY OF COMBAT 

F. W. Lanchester (1868-1946) did pioneering work in 
military operations research by attempting to quantitatively 
justify VonClausewitz ' s principle of concentration. His 
work led him to postulate a simple deterministic model for 
the mutual attrition of two homogeneous forces. This 
pioneering work first appeared in a series of articles in 
British Journal Engineering during 1914. Since that time 
deterministic formulations have been vastly extended. The 
attrition process has also been viewed from the stochastic 
standpoint. 

During World War II, B. Q. Koopman extended Lanchester 's 
original results and also suggested a reformulation of the 
problem in stochastic form. This stochastic attrition 
process has been approximately termed by Koopman as the 
Lanchester stochastic process. Other workers have subse- 
quently employed a stochastic analysis of combat. 

Among these various works, Brackney in his paper, "The 
Dynamics of Military Combat," ([3]), categorized the types 
of combat situation into nine separate, clearly defined 
situations . 

Subsequent workers (including Smith [8] whose work is 
reviewed below) have used Brackney 's scheme for classifying 
combat situations. 



5 



Analogous to the deterministic solution of Lanchester's 
equation, Koopman (whose work is reported in Morse and 
Kimball [5]) developed the probability that there will be 
some survivors on one side when the other side is anihilated 
(i.e. the probability of winning) when the casualty rate on 
each side is proportional to the product of the number of 
firers and the number of targets. This, of course, cor- 
responds to the deterministic "linear-law" attrition process. 

Later, Brown [6] developed the probability of winning 
corresponding to a deterministic "square-law" attrition 
process (i.e. casualty rate on each side proportional to 
the number of enemy firers) . Brown also gives some approx- 
imations for computing the probability of winning. Smith 
[8] gave the solution for the remaining case ("mixed law"). 

In this thesis results from the stochastic and deter- 
ministic Lanchester models are contrasted, much in the same 
spirit as George H. Weiss did in [9]. 

B. STUDY OBJECTIVES 

The term "Lanchester theory of combat" is commonly taken 
to refer to a set of theories which attempt to explain 
changes in force levels due to combat attrition in terms of 
weapon system performance characteristics, force levels, 
and composition. There are two general approaches to 
modeling the combat attrition process: 

(1) deterministic formulation which takes the form of 
a system of first order differential equations 

(2) Stochastic formulation v/hich viev;s the casyaktt 
process as a Markov process. 
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In this thesis some idealized models of combat between 
two homogeneous forces are considered. The forecasts of 
victory from stochastic and deterministic models are com- 
pared. Numerical computations of the probability of winning 
are performed according to theoretical results which have 
appeared in the literature. 

Numerical computations of probability of winning were 
done by the IBM 360/67 computer at the Naval Postgraduate 
School. The computational results will be discussed by 
considering graphs and tables for various cases. 

Similarities and differences between the deterministic 
and the stochastic models will be discussed. Based upon the 
computational results obtained, the tactic "Don't attack 
an enemy position without at least a three-to-one initial 
force ratio" will be evaluated 

The model parameters selected for the computations 
presented in this thesis were chosen to be in agreem.ent 
with the military experience of the authors. Although 
rather idealized situations were studied in order to preserve 
analytic tractability of the models, many of the quantita- 
tive results discussed in this thesis apply in more complex 
situations . 
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II. COMBAT BETWEEN TWO HOMOGENEOUS FORCES 



Brackney [3] considered nine types of combat situations. 
These nine situations correspond to all the ways of combin- 
ing three basic attrition-rate schemes for each side. Three 
basic situations are shown below for combat between homo- 
geneous X- and Y-forces. 



dx 

dt 



- a 




dt 



-ay ; 




dx 

dt 



-axy ; 




( 1 ) 

( 2 ) 

(3) 



where x and y are the numbers of surviving forces at time 
t on Red and Blue side respectively, a and 3 are the constant 
of proportionality (attrition-rate coefficients) and are 
equal to the effectiveness per unit force for each side. 

We shall call (1) the linear law, (2) the square law, 

(3) the mixed law. The reason for the name "linear," 

"square" and "mixed" law is explained by Lanchester [2]. 

For mathematical convenience, only the case of fighting 
to finish, that is when one side is anihilated , is con- 
sidered in this thesis. 



A. LINEAR LAW ATTRITION ON BOTH SIDES 
1 . Assumptions 

(1) Two forces attack each other. Each unit (man) 
on each side is within range of all units of the other side. 
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(2) Units on each side are identical but the units 
on one side may have a different attrition rate than the 
opposing units. 

(3) Each firing unit is aware only of the general 
area in which enemy forces are located and fires into this 
area without knowledge of consequence of this fire. 

(4) Fire from surviving unit is distributed uni- 
formly over the area in which enemy are located. 

2 . Deterministic Model 

The two sides are designated as "Red" and "Blue". 

X , y = initial force levels. 

>^(t), y(t) = force levels on the side of "Red" and 
"Blue" respectively at time t. 

Then the Lanchester equations are: 



dx 

dt 



-axy 




-3xy 



The solution with time eliminated is 
3{x^-x(t)} = a{y^-y(t)} 

Since we are considering the fight to finish only, 
we could think the probability of winning for the Blue can 
be represented as follows: 



P(B,y^) 



X 

1 if a/e > — 

— Y 
^ o 

0 if a/e <_ ^ 

^o 



where P(E,y^) is the probability of Blue v'in with initial 
force level y^. This probability can be charged as 
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if 1 > 



P(B,y^) = 



ay 



3x, 



if 1 < 



— ay 



o 
3x 

Let us hence forth refer to as the normalized 

ay 



force ratio. 

3 . Stochastic Model 

Define P(0,n,M,N) as the probability that there will 
be n survivors on the Blue side when the Red side is 
anihilated . 



P(B,N) the probability of Blue win with initial 
force level N 



where M,N = initial numbers of units 

m,n = number of surviving units on the side of "Red" 
and "Blue" respectively at time t. 



Then 

P(0,n,M,N) = (^) P(0,n,M-l,N) + (^) P(0,n,M,N-l) 

ct ot 

for notational conveneince, let p = — r-pr > <3 = — rn a-nd 

a+3 cx+3 

P(0,n,M,N) = F(M,N) 

F(M,N) = pF(M-l,N) + qF(n,N-l) 

with F(0,n) = 1 
F(0,n) = 0 



otherwise solution (can be verified by substitution) is 

/M+N-n-1 



F(M,N) = 



M-1 



M N-n r, . 

p q = P (0,n,M,N) 



10 



Hence 



N 

Z P(0,n,M,N) 
n=l 



P(E,N) = 



N /M+N-n-l\ „ ,, 

= E p q 

n=l \ M-1 / 



let j = N-n 

then n = 1 implies j = N-1 
n = N implies j = 0 



Then 



N-1 /M+j-1' 
P(B,N) = Z 

j=0 \ M-1 



M j 
p q-' 



or 



P(P,N) = 



/ __?L 



n 



N-1 /M+j-1' 



A I : 



a+3 



B. SQUARE LAW ATTRITION OR BOTH SIDES 

1. Assumptions 

(1) as in Case A 

(2) as in Case A 

(3) Each firing unit is sufficiently well aware of 
the location and condition of all enemy units so that when 
a target is killed, fire may be immediately shifted to a 
new target 

(4) Fire is uniformly distributed over surviving 

units , 

2 . Deterministic Model 

The Lanchester-type equations are 
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dt 



= -ay 



^ = -6x 
dt 



with solution 



e{x^-x^(t)} = a{y^-y^(t)} 



Then the probability of winning for Blue can be 



represented as: 



P(B,y^) = 



if a/B ^ 

yo 



P(B,y^) = 



0 if 


a/B 


^o 




changed as 


1 if 


1 > 


/B/a 


^o 

(-^) 

Yo 


0 if 


1 < 


/B/a 


Xn 

{—) 

Yo 



We refer to /B/a ( — ) as the normalized force ratio 

^o 

6 ^o"* 

and — — =“ as the squared normalized force ratio for square- 

Ot 

i'o 

law attrition. 

3. Stochastic Model 



P(0,n,M,N) = (-S^)P(0,n,M-l,N) + P (0,n,M,N-l) 



'aN+BM. 

let F(M,N) = P(0,n,M,N) 

Then solving for 



aN+BM 



’ (M,N) 
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with 



F(0,n) = 1 
F(0,N) = 0 



otherwise we can get [6] 




and 




P(B,N) = 2 ^ 

k=0 (N-k) ! k! r (M + — + 1) 



C. THE MIXED CASE 
1 . Assumption 

(1) As in Case A 

(2) As in Case B 

(3) As in Case A for Blue side and Case B for Red 



side 



2. Deterministic Model 



The Lanchester-type equations are: 




with solution 



2g {x^ - x^(t)} = - y^(t)> 



Then the probability of winning for Blue can be 



represented as: 




0 



1 



if 



if 




Yo 
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This can be changed as 



Define 

3. 



P(B,y^) 



1 


if 


1 > 


0 


if 


1 < - 



2$x 
c 

2 

Vo 

2 3Xq 



23x 



as normalized force ratio for mixed case. 



Yo 

Stochastic Model 

P(0,n,M,N,)= (0,n,M-l,N) + ( P (0 , n , M, N-1) 



after solving this 



P (0,n,M,N) 



N 

E 

£=n 






M 



^N-n-1 . ^ 
r(N- +1) r (n-n+i) 



and 



P(B,N) 



11=0 U+-)*’ r (N-t+i) r (1) 

a 
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III. COMPUTATIONAL RESULTS 



A. LINEAR LAW 

A digital computer program was developed to perform the 
rather tedious computation of the probability of winning 
for various levels of model parameters. 

Plots of probability of winning for Blue versus nor- 
malized force ratio are given on pages 19-24. 

In order to clarify the probability shift, each graph 
in Figures 1 to 5 are drawn for different initial force 
level of Blue forces. In Figure 6 the quantity M+N is held 
constant. In Figures 1 to 5 the results of using different 
values of the attrition-rate coefficient can be seen. The 
straight line of each graph represents the probability of 
winning for deterministic la\^. 

As we can see in these figures, increasing the number of 

combatants causes the probability of winning as computed 

from the stochastic model to approach the deterministic win 

criterion (i.e. x wins when Bx > ay ) . 

o o 

It is reasonable to use the deterministic model to study 
combat betv/een two forces v/ith large numbers of combatants. 
This is because the mathematical convenience obtained from 
using the deterministic model is much greater than the pos- 
sible errors between the two models. 

One really interesting result is that conditions to win 
for both deterministic and stochastic models turn out to be 
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very similar. In case of deterministic model if the nor- 
malized force ratio is less than 1, the probability of 
winning for Blue is 1 all the time. This implies Blue side is 
always sure to win. But in case of stochastic model, when 
the normalized force ratio is less than 1, there is still a 
possibility to loose the battle. Especially when the number 
of Blue combatants are less than 40, the difference between 
deterministic is critical. An extreme case is the example 
from Figure 1 v/hen N = 2 and M = 1. In this case we can see 
the probability of winning for Blue is equal to 1 for deter- 
ministic case. But from the stochastic model the winning 
chance is only 3/4. This means the Red side still has 1/4 
chance of winning. 

When the initial number of Blue combatants is less than 
40, the critical part for the normalized force ratio is 
between 0.5 and 1.5. This implies that for the combat, in 
which the relatively small number of combatants are partici- 
pating and the normalized forces ratios are . 5 apart from 
unity, using the stochastic model is reasonable. But for the 
case when the large number of combatants are participating 
or normalized force ratio are large enough, using the determinis- 
tic m.odel is m.ore realistic. 

Another interesting result is a probability shift 
apparently caused the attrition-rate coefficient difference. 

In Figures 2 to 5 when the normalized force ratio is unity 
the probability of winning is different from .5. Before 
performing these computations it was anticipated that the 
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probability of winning would be .5. But these numerical 
results showed that this guess was not true. Also, the 
probability difference from 0.5 becomes smaller as the 
initial number of combatants becomes larger. The attrition- 
rate coefficients influence the probability of winning more 
than increasing the initial number of combatants. 

Figure 6 also indicates when larger number of forces are 
involved, the probability that the larger force will win 
approach unity. The probabilistic form is preferable to the 
simple differential equation, however it is not clear that 
the additional difficulty of solution is consistent with the 
improved realism in view of the other known variable of 
actual combat. 

It seems appropriate to examine the often quoted maxim 
of military tactics that one should not attack without at 
least a three-to-one force ratio. From Figures 1 to 5 , the 
region where the normalized force ratio is less than 1/3, 
the probability of winning is greater than .99. This implies 
Blue side is more than 99% sure to win. But it is natural 
result because the Blue side has 3 times greater "combat 
power" than the Red side. 

Hence we may investigate another v/ay when 

aN > BM Blue win 

if N = 3M the winning condition become 
3aM ^ BM 

This can be written 
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To analyze this condition we made Table I. From column 
6 of Table I we can see the probability of winning is 
inversely proportional to force ratio and approach to 0.5 
when the number of combatants for Blue side are increasing. 
This implies in linear law battle if ^ = 3, then the 3 to 1 
force ratio attack gives the Blue side at least draw. 

Column 8 shows the winning chance of the Blue side is very 

g 

small because the winning condition — _< 3 is not satisfied. 

B . SQUARE LAW 

The same procedure for the analysis using "linear-law" 
attrition has been done for square-law attrition except the 
abcissa of the graph of square law is the square of nor- 
malized force ratio. 

The computational results are very similar to linear 
law case. To avoid the same citation twice, we try to 
analyze the square law case by comparing the difference 
from linear law case. 

One difference is the probability shift which cause from 
the different value of the ratio of attrition rate. In the 
case of square law the attrition-rate coefficient influence 
less the probability of winning than in the case of linear- 
law attrition. 

To investigate the most usual situation in combat, where 
one antagonist attacks another; v;e can develop condition of 
winning by considering the practical forces ratios. 
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Figure 1 



^ Linear Law: a = 0.05, 6 = 0.05 
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Linear Lav;: 
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ct = 0.2, 3 = 0.5 
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Figure 3 
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Figure 4 

Linear Law: a = 0.3, '3 = 0.1 



P(Y . ) 

win 
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Figure 5 

P(Y . ) 

win Linear Law; a =0.1, 8=0.3 




23 



tu- 



^ .0 









n 



I. 



i 



Figure 6 



P(Y ■ ) Linear Law: a = 0.5, B = 0.5 
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2 2 
aN ^ 6M 

if N = 3M the winning condition will be 
a ( 3M) ^ 

This can be written 

■^ < 9 
a — 

By investigating the column 4, 1, 9 of Table I, we 
can expect the same result as linear law case. 

C. MIXED LAW 

All results are similar to those for cases A and B above. 
From column 5 of Table I the probability of winning is small 
compared to case A and B. This is caused from the effect of 
idealized assumption for mixed law. The actual situation 
of combat is quite different from the idealized situation, 
attacker uses area fire while the defender use aimed fire. 

Below the line where the initial force ratio 8 to 24, 
the probability of v/inning become much greater. This shows 
the difference between the deterministic and stochastic 
model more clearly. 
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Figure 7 
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Figure 8 
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Figure 9 



Square Lav;: a = 0.2, 6 = 0.05 
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Figure 10 
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Figure 11 



P (Y ) 

win Square Law: a = 0.1, 3 = 0.3 
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Figure 12 
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Figure 13 



P (Y . ) Mixed Law: a = 0.05, R = 0.05 

Win ^ 
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IV. SUGGESTED FUTURE WORK 



The stochastic model solution includes many factorial 
and gamma functions. This apparently gives computational 
difficulties due to the large numbers generated causing 
"overflows" in the computer. To overcome such difficulties 
recursive formulas were developed for the case of linear- 
law model. These are shown in the following table. 



TABLE II 



N 


P (B,N) 


1 


d” 


• 1 


2 


D”j 


1+MC j 


3 


D«! 


[l+MC + C^J 


4 


d” 


C2 + M(M+1)(H+2^ c 2 j 


• 




• • 

# • 



a+6 ' ^ a+B 

Because the recursive formulas for the case of square 
and mixed model were not developed at this time, another 
approach had to be taken to eliminate overflow problems for 
the "Square-Law" and "m.ixed-law" models. This v/as 
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accomplished by using logarithms in the appropriate formulas 
for the probability of winning. Although this took much 
computing time and was very inefficient, it did allow results 
to be computed for force levels up to approximately 20 to 50. 
Development of such recursive formulas for square law and 
mixed law cases was left for future work. 

Secondly the current work examined only the simplest 
cases of combat between two forces in the Lanchester theory 
of combat. If we introduce the time to the model then the 
situation becomes more complicated. We can expect more 
interesting results by introducing reality into the model 
and do the same work. 
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V. CONCLUSION 



Idealized combat situations in the Lanchester theory 
of combat were examined using both deterministic and sto- 
chastic models. In the simplest case of combat between two 
homogeneous forces, the probability of winning was studied 
and results contrasted for these two types of models. Based 
on the results of numerical computation of the probability 
of winning, it has been concluded that the deterministic 
model (even though win probabilities are either zero or one 
depending upon whether the "deterministic win criterion" 
is satisfied) yields satisfactory forecasts of comJDat out- 
comes (probability of winning) when there are large numbers 
of combatants (at least on one side) . The curve of the 
probability of winning (in a fight-to-the finish) versus 
the "normalized force ratio" (or appropriate equivalent) 
is characterized by a steep slope in the probability of 
winning for normalized forces ratios between 0.5 and 1.5. 
This slope becomes steeper as the number of combatants 
increases. Thus, in such circumstances the addition of a 
relatively few additional forces on one side may signifi- 
cantly increase that side's chances of winning. 

The results of tedious computation shows the attrition- 
rate coefficients influence more the probability of winning 
than increasing the initial number of combatants. The three 
to one force superiority is needed to win the combat with 
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much confidence. The results of three cases 
square law and mixed law, are very similar, 
is much work to be done in the future. 



i 



linear law, 
Still there 
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COMPUTER PROGRAM 



COMPUTER PROGRAM FOR LINEAR-LAW 



IMPLICIT REALM'S! A-H,0-Z,$) 

9000 t^ORMAT (* » , 5X , F9 . 5 , 5X , F9 . 5 , 5X , F9 . 5 , 5X , F o . 5 , 5 X , F 9 . 5 , 

^5X,I5,5X,I5) 

9100 FORMAT (» • , 8X , • A » , HX , • B ' , lOX , • A/B ' ,9X , ' 3M/AN' , 
*5X,’PR0B WIN' , 7X, 'M',7X, 'N • ) 

WRITE (6,9130) 

A=0.05 
B=0.05 
AOE=A/R 
C=B/( A+B ) 

D=A/( A + B ) 

DO 1300 N=1,133 
DO 1000 M= 1,130 
XM=nPLOAT (M ) 

XM=DRLOAT(N) 

XM0N=(B*XM )/( A*XM) 

TERM1=0**M 

SUM=1.0 

TFRM=1.0 

IF(N.LT.2) GO TO 200 
DO 100 1=2,N 
XI=DFLCAT ( I ) 

M2I=M-2 + I 

XM2I =DFL0AT{ •■<21 ) 

TFRM = TERM*XM2 1/ ( XI-1.0 )="C 
SUfM= SU'-' + TFRM 
100 CONTINUE 
200 PR0B7=TERM1*S JM 

WP I TE ( 6 , 90 3 3 ) A , B , AOB , X NON , P ROPY , M , N 

1000 CONTINUE 
STOP 
END 



COMPUTER PROGRAM FOR SQUARE-LAW 



IMPLICIT REAL’^3 ( A-H,0-Z,$ ) 

9000 FORMAT(' • , 5X , “ 9 . 5 , 5X , f 9. 5 , 5X , R R . 5 , 5 X , ^ 9 . 5 , 
<^5X, F9.5 ,5X, IE, EX, 15 ) 

A = 0.3 
B = 0. 1 
A03=A/B 

DO 1000 N=2,130 
XN=DFL0AT (N ) 

DO 1000 M=2,130 

XM=3FL3AT(M ) 

MN=M+N 

XMN=DFLOAT (MN) 

XMON= ( B*XM^*2 ) / { A^X^'^'^Z) 

I F (XMON.LT .0 .2 ) GO 100 3 
IF( XMQM.3T. 13. 3) 30 TO 1333 
TFRM1=XM«DL0G( AOB ) 

SUM=0.0 
00 100 K=1,N 
XK=DFLOAT( K ) 

TERM2=3. 0 
TERM3=0. 0 
TERM4=3 .0 
NHI NK=N-K 
DO 111 K1=1,NMINK 
XK1=mfL 0AT (Kl) 

TFRM2=TERM2 + '^L0G{XK1 ) 

111 CONTINUE 

DO 112 K2 = l,’' 

XK2=DFL0AT( K2) 
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TFRM3 = tFRM3+CLD3 ( XK.2 ) 

112 COMTIMUF 
XKAB = XK^-AOF 
DO 113 <3=1, M 
XK3=0FL0AT (K3 ) 

XK3XK=XK3+XKA5 

T E P. M 4 = T E P M A + 0 L 0 G { X K 3 X K ) 

113 CONTINJ= 

XJC0E = (-1.0)’"*=MMIMK 
XKM0S = X^»M*DL03 { XK ) 

STUB=TERM1 + XK.MMS-TER02-TERM3-T ERMA 
SUM=SU^1 + XJC0E<=DEXP{ STUB) 

100 COMTIMJE 
PR0BY=SUM 

WR! TE( 6, 9000) A,B,AD? , XMOM ,PRCBY , M, N 
1000 COMTIMJE 
STOP 
END 



COMPUTER PROGRAM FOR MIXED-LAW 



IMOLIGIT REAL*^8( A-H,0-Z,$) 

9000 FORMAT!' ’ , 5X , P 9 . 5 , 5X , F9 . 5 , 5 X, P 9 . 5 , 5X , 9 . 5 , 
«5X,F9. 5, 5X, 15, 5X, 15 ) 

A=0. 05 
B=0.05 
A0B=A/B 
BOA=B/^ 

DO 1000 0=2,100 
XN = DPL0AT( NO 
DO 1000 M=2,100 
XM=DFLCAT (M) 

MN = '^+-N-1 

XM0N=l 2. O'-B-XM )/ ( A*XN='-=2) 

XMN=PFLOAT (MK) 

SUM=0. 0 
DO 100 L=1,M 
XL=D'=L0AT{ L) 

term 1 = 0. 0 

TERM2 = 0 .0 
NML=M-L 

XNML = DFL':'AT( MML) 

LM1=L-1 

XLMl =DFLOAT ( LMl ) 

XJC0‘= = { - 1. 0) *^NML 
IFINML.LP.l) 30 TO 200 
DO 111 L1=1,NML 
XL1=DFL0AT(L1) 

TFR»^1=TERM1+DL0G(XL1 ) 

111 CONTINUE 
GO TO 300 

200 T=RM1=0.0 

300 lE(LMl.LE.l) 30 TO AOO 
DO 112 L2=1,LM1 
XL2=DFL0AT {L2 ) 

T‘=RM2=TFPM2 + DL03 (XL2 ) 

112 COMTIMJE 
GO TO 500 

AOO T^RM2=0.0 
500 XLBCA=XL+8CA 

STU5 = XMM*0LCG( XL )-XM=DLDG(XLBOA )-T ERMl-T ERM2 
SUM=SU'^ + XJCOP<'OEXP( STUB) 

100 CONTINUE 
PROBY=SlJM 

WRITP(5,9000) A,B,A03 , XMOM , PR OB Y , M , N 

1000 COMTIMJE 
STOP 
END 
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